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C$ • 1 . Introduction vide a direct way to trace the bulk of the star formation in the 

early universe. 

Pointed submillimeter and millimeter observations of opti- 
cally selected quasars or radio galaxies have the advantage over 
I20031 iBertoldi et ai1l2003al) have provided a view of the dust deep field surveys in providing unique source identification and 
content of a few hundred galaxies and quasars in the redshift redshifts. In particular, targeted observations of optically lumi- 
range 1 < z < 6.4. In these objects, a significant fraction of nous quasars have been successful in revealing massive bursts 
the energy generated by star formation is processed by dust of star formation in their host galaxies. Although optically lu- 
and re-emitted at far-infrared wavelengths. The average space minous quasars, which trace the most massive collapsed struc- 
density of ultraluminous infrared galaxies and quasars at high-z tures to have formed in the early universe, are rare objects 
is found to be thousandfold greater than in the local universe, which bear little relation to the objects that make up most of the 
Observations at millimeter and submillimeter wavelengths pro- far-infrared and submillimeter background, they provide pow- 
erful probes to study the relation between star formation, mas- 
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Abstract. We report the detection of CO emission from the optically luminous, radio-quiet quasar J140955. 5+562827 (here- 
after J1409+5628), at a redshift z C o = 2.583. We also present VLA continuum maps and VLBA high spatial resolution observa- 
tions at 1.4 GHz. Both the CO(3— >2) and CO(7— >6) emission lines are detected using the IRAM Plateau de Bure interferometer. 
The 3— »2/7— »6 line luminosity ratio is about 1/3, indicating the presence of warm and dense molecular gas with an estimated 
mass of 6x 10 10 M G . The infrared-to-CO luminosity ratio Z-fir/£'cq(i->0) ~ 500 Lo (K km s -1 pc 2 )~', comparable to values found 
for other high-z sources where CO line emission is seen. J1409+5628 is detected using the VLA with a 1.4 GHz rest-frame 
luminosity density of 4.0 x 10 25 WHz -1 . The rest-frame radio to far-infrared ratio, q, has a value of 2.0 which is similar to the 
values found in star forming galaxies. At the 30 mas resolution of the VLBA, J1409+5628 is not detected with a 4<r upper limit 
to the surface brightness of 0.29 mly beam -1 . This implies a limit to the intrinsic brightness temperature of 2 x 10 5 K at 8 GHz, 
typical for nuclear starbursts and more than two orders of magnitude fainter than typical radio-loud active galactic nuclei. Both 
the properties of the CO line emission and the radio emission from 11409+5628 are therefore consistent with those expected 
for a star forming galaxy. In 11409+5628 young massive stars are the dominant source of dust heating, accounting for most 
of the infrared luminosity. The massive reservoir of molecular gas can sustain the star formation rate of a few 1000 M Q yr~' 
implied by the far-infrared luminosity for about 10 million years. 

Key words, galaxies: formation - galaxies: starburst - galaxies: high-redshift - quasars: emission lines - quasars: individual: 
1140955.5+562827 - cosmology: observations 



In the recent decade, millimeter 


and submillimeter deep blank 


field surveys (e.g., 


Ivison et al. 


2000; 


Blain et al. 


2002) and 
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* This paper is based on observations obtained with the IRAM |20Q3bl[Walter et all20 03) 
Plateau de Bure Interferometer. IRAM is funded by Centre National 

de la Recherche Scientifique (France), the Max-Planck-Gesellschaft The study of the redshift range 2 < z < 3 is important 

(Germany), and the Instituto Geografico Nacional (Spain). in studying the star formation history in the universe since it 
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traces th e peak of the space density of the quasar population 
JShaver et alll996l) . an d corresponds to the m edian redshift of 
submillimeter galaxies ( Chapman et alJ 120031) . Recent pointed 
observations of the 1.2 and 0.85 mm thermal dust continuum 
emission of z ~ 2 optically luminous, radio-quiet quasars re- 
veal that 1/3 of these quasars are also luminous in the infrared 
with far-infrared luminosities of Lfir ~ 10 13 L and esti- 
mated dust masses of typically ~ 1O 8 M dOmont et al.ll2003t 
IPriddev et alJ2003l) . 

The optically very bright {Mb = -28.4), radio-quiet 
quasar J1409+5628 is by far the strongest mm source in the 
Omo nt et alJ l2003l> survey. With an estimated far-infrared lu- 
minosity of 4 x 10 13 L , J1409+5628 ranges amongst the 
most luminous infrared high-z sources found to date. Under the 
assumption that the dust is predominantly heated by massive 
stars, the inferred s tar formation rate is several 1000M o yr _1 
( Om ont et al 12003). However, the fraction of infrared luminos- 
ity due to dust heated by star formation and dust directly heated 
by the Active Galactic Nuclei (AGN) in J 1409+5628 is still an 
open issue, as is the case for almost all of the high-z quasars 
studied to date. Both the large masses of dust and the high star 
formation rate indicate that J 1409+5628 should have a copi- 
ous reservoir of molecular gas to sustain the formation of a few 
1000 M per year over its dynamical time. 

In high-z sources, the search for molecular gas is best done 
by using the redshifted rotational lines of CO. The molecular 
gas in star-forming galaxies is warm and dense enough to excite 
the higher CO rotational levels (J > 3) which are shifted into 
the millimeter windows for redshifts greater than 2. However, 
a good determination of the redshift of the target object, either 
from optical or near-infrared spectroscopy, is a prerequisite for 
a detection because the bandwidths of the current heterodyne 
receivers are still relatively narrow. For example, at the Plateau 
de Bure Interferometer, the available instantaneous bandwitdh 
is Av = 580 MHz which corresponds at 96 GHz to a velocity 
range of 1810 km s _1 or Az ~ 0.02 at z = 2.6. 

Despite this current limitation, over the last few years CO 
line emission has been detected in 24 far -infrared luminous 
high-z quasars, radio and (sub)mm galaxies (Ca rilli et alJ2 004. 
and references therein). These observations reveal large reser- 
voirs of molecular gas (10 11 M ), a prerequisite for efficient 
star formation, and prompted also the first studies on the prop- 
erties of the dense interstell ar gas in starburst galaxies in the 
early universe (e.^. | Barvainis et all 1 9971: ISolomon et all2003l 
iBertoldi et alJl2003bl) . In addition, these observations allow us 
to estimate, unhindered by extinction, the dynamical masses 
in these systems, whic h are the k ey constraining hierarchical 
models (see, e.g. lNeri et all2 003i. 

In this paper, we present a study of the CO emission in 
J 1409+5628 made with the Plateau de Bure interferome- 
ter together with observations of the radio emission carried 
out with the Very Large Array and the Very Large Baseline 
Array. The data provide strong evidence that the origin of 
the infrared luminosity is dominated by the starburst activ- 
ity of the host galaxy of the QSO J1409+5628. Throughout 
this paper, we adopt the concordance A-cosmolo gy with Ho = 
71 km s- 1 Mpc" 1 , Q A = 0.73 and Q m = 0.27 - ispergel et alJ 



120031) . At z = 2.58, the luminosity distance Di is 21.5 Gpc and 
an angular scale of 1" corresponds to 8.1 kpc. 

2. Observations 

2.1. Plateau de Bure Interferometer 

Observations were made with the IRAM Plateau de Bure inter- 
ferometer in May and June 2002 and from January to March 
2003. We used the D and B configurations, as well as non- 
standard configurations, with 4 to 6 antennas, resulting in a 
good coverage of the uv plane with baselines ranging from 
24 to 400 meters. This results in final synthesized beams 
of 2" 44 x l'/65 at a position angle of 60° at 3.1 mm, and 
1.01" x 0.55 " with a position angle of 52° at 1.3 mm. 

The dual frequency set-up was used to search simultane- 
ously for the CO(3— >2) emission line and dust emission at 3 
and 1.3 mm, and subsequently for the CO(7— >6) emission line. 
During the first series of observations, the 3.1 mm receiver was 
tuned to a central frequency of 96.780 GHz, corresponding to 
the CO(3->2) line at a redshift of z = 2.5730, close to the 
strong Lyman-a absorption line seen at z = 2.5758 in the op- 
tical spectrum of J1409+5628 (see Sect. ED. The CO(3^2) 
line emission was tentatively detected at the low frequency end 
of the available bandwidth. In subsequent observations, we re- 
tuned the 3.1 mm receivers to 96.505 GHz, corresponding to a 
redshift of z = 2.583 and confirmed this tentative detection. 
After the first tentative detection of the CO(3— >2) emission 
line, we retuned the 1.3 mm receivers in the upper sideband to 
225.198 GHz corresponding to the redshifted CO(7— >6) emis- 
sion line, while the lower sideband was tuned to a frequency 
3 GHz away (at 225 GHz) to measure the 1.3 mm continuum. 
Typical SSB system temperatures were « 150 K and » 400 K 
at 3 mm and 1.3 mm respectively. 

The data were reduced, calibrated and analyzed using the 
IRAM package Gildas. During the data reduction, we flagged 
antennas with too low efficiencies or too high system temper- 
atures, and discarded data taken during bad atmospheric con- 
ditions. The final usable datasets correspond to total on-source 
integration times of 52 and 27 h at 3. 1 and 1 .3 mm, respectively. 
After phase and amplitude calibration of each session, we com- 
bined the data from the different settings into one single table 
at full spectral resolution. The flux calibration was done using 
standard calibrators 3C 273, CRL 618 or MWC 349. 

2.2. Very Large Array 

J1409+5628 was observed with the Very Large Array (VLA) 
of the NRAO 1 at 1.4 GHz in the A configuration (maximum 
baseline of 30 km), with a total bandwidth of 100 MHz and two 
polarizations. The total observing time was 2 hours. The en- 
tire dataset was self-calibrated using field sources, and standard 
phase and amplitude calibration was applied for J1409+5628. 
The absolute flux density scale was set with the observations of 
3C 286. 

1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement 
by Associated Universities, Inc. 
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Fig. 1. Left panel: Velocity-integrated map of the CO(3— »2) line emission toward J 1409+5628. The contour step is 0.3 mJy/beam, 
corresponding to 2cr, dashed contours are negative. Right panel : Velocity-integrated map of the CO(7— >6) line and 1.3 mm 
continuum emission toward J 1409+5628. The contour step is 2 mJy/beam, corresponding to 2cr, dashed contours are negative. 
In both panels, the cross shows the optical position of the quasar taken from the 2MASS catalogue. The synthesized beams are 
shown in the lower left corner of each panel. 



The final image was generated using the wide field imaging 
and deconvolution capabilities of the NRAO's Astronomical 
Image Processing System (AIPS) task IMAGR. The present 
1.4 GHz continuum data achieved a rms noise (<x) of 16/iJy 
which corresponds to the theoretical expected noise level. The 
Gaussian restoring CLEAN beam full width at half-maximum 
(FWHM) is ~ 175. A detailed discussion of the VLA data will 
be given in Petric et al. (in preparation). 

2.3. Very Long Baseline Array 

We observed J1409+5628 with the VLBA of the NRAO on the 
9th of March and 7th of April, 2003, for a total of 14 hours us- 
ing the standard continuum mode at 1 .4 GHz (total bandwidth 
of 16 MHz in two dual circular polarizations). The calibrator 
J1408+5613 was used for phase referencing with a calibration 
cycle time of 3.5 minutes. Standard a priori gain calibration 
was applied. At 1 .4 GHz, the VLBA short spacing limit filters 
out all spatial structure larger than about 0715. 

A number of test cycles were also included to moni- 
tor the coherence of the phase referencing. These tests in- 
volved switching between two calibrators (J1408+5613 and 
J1419+5423) with a similar angular separation and cycle time 
as that used for the target source. Images of the second 
calibrator (J1419+5423) were deconvolved using two differ- 
ent approaches: by applying the phase and amplitude self- 
calibration solutions of the phase reference source J1408+5613 
on J1419+5423, and by self calibrating J1419+5423 itself, both 
in phase and amplitude. The ratio of the peak surface bright- 
ness between the final images of the two approaches gives a 
measure of the effect of residual phase errors after phase ref- 



erencing (i.e. 'the coherence' due to phase referencing). At all 
times the coherence was found to be better than 95%. Also, 
phase referencing as used herein is kno wn to preserve ab solute 
astrometric positions to better than 071 (Fomalo ntl 19991) . 

3. Results 

3.1. Molecular gas 

The velocity-integrated maps of the CO(3— >2) and CO(7— »6) 
emission lines of J1409+5628 are shown in Fig. and the 
corresponding spectra are displayed in Fig. [2 The CO(3— >2) 
line is detected with high signal-to-noise ratio at a central fre- 
quency of 96.504 GHz, corresponding to a redshift zco = 
2.5832 ± 0.0001. At the same redshift, the CO(7^6) emission 
line is also detected on top of the 1.3 mm continuum. 

Within the astrometric uncertainties, the CO emission and 
the 1.3 mm continuum, centered at 14 h 09 m 55 s :5 + 56°28'2674 
(J2000.0), coincid e with the optical pos ition listed in the 
2MASS catalogue (Bar khouse & Halll200lh . At the 1.0"x0.5" 
resolution of our Plateau de Bure observations at 1.3 mm, the 
CO(7— >6) emission in J1409+5628 is still unresolved, yield- 
ing an upper limit to the linear size of the region emitting in 
CO of5kpc. 

We fitted the high signal-to-noise CO(3— »2) spectrum with 
a gaussian profile where the central frequency (v bsX the full 
width half maximum (AvfwhmX the continuum flux density 
(S v ), and the integrated CO emission flux density (Ico) were 
kept as free parameters. The peak intensity was derived from 
the gaussian fit parameters. The best gaussian fit yields a line 
width to the CO(3— >2) emission of 311 ±28 km s " 1 , sim ilar 
to the widths found in other high-z quasars (e.g., ICox et all 
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Table 1. Properties of the CO lines observed toward J 1409+5628 



Line 




Zco peak int. 


Avfwhm 
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^co 


L co 


Leo 




[GHz] 


[mJy] 


[kms" 1 ] 


[mJy] 


[Jy km s -1 ] 


[10 10 Kkms- 1 pc 2 ] 


[10 8 L ] 


CO(3^2) 


96.504 


2.5832 ± 0.0001 6 ± 1 


311 ±28 


< 0.5 1 


2.3 ± 0.2 


8.2 + 0.6 


1.1 ±0.1 


CO(7^6) 


225.1197 


10 ±3 




6±2 


4.1 ± 1.0 


2.6 ±0.7 


4.4 ± 1.2 



Notes - f 3o" level upper limit. Adapting the CO(3— >2) parameters (see text) 
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Fig. 2. Observed spectra of the CO(3— >2) and CO(7— >6) line emission toward J1409+5628, with velocity resolutions of 58 and 
76 km s _1 , respectively. The offset velocity corresponds to central frequencies of 96.504 and 225.1197 GHz. The error bars 
represent the statistical uncertainties from the fit in the Mv-plane. The dashed lines show the gaussian fits together with the fit 
to the continuum with the parameters listed in Tabled The dot-dashed line is the fit to the CO(7— >6) line emission leaving the 
linewidth as a free parameter 



2002), and a peak intensity of 6 + 1 mJy. The velocity- 
integrated CO(3— >2) flux is 2.3 + 0.2 Jy km s _1 . At the position 
of J 1409+5628 we obtain from the combined fit a 3<x upper 
limit to the continuum flux density at 3.1 mm of Sye.sGHz < 
0.5 mJy. 

The width and centroid of the gaussian fit to the CO(3— >2) 
line were adopted for the profile of the lower signal-to- 
noise ratio CO(7— >6) emission line. The peak intensity of 



the CO(7— »6) emission is 10 + 3 mJy, the velocity-integrated 
CO(7— >6) flux 4.1 ± 1.0 Jy km s~', and the 1.3 mm continuum is 
detected at a flux density of 6 ± 2 mJy. The CO(7— »6) emission 
line is detected with a signal-to-noise ratio of 4, as compared 
to the ratio of 1 1 in the case of CO(3->2). The CO(7^>6) line 
flux is therefore somewhat uncertain, and the precise value de- 
pends on how the spectrum is analyzed. For instance, leaving 
the line width as a free parameter in the gaussian fitting yields 
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Avfwhm = 229 + 84 km s -1 and a CO(7— >6) integrated line flux 
of 3.0 ± 1.3 Jy kms -1 , i.e. a 2.3 cr result (Fig.|2jl. 

The upper limit of the continuum flux density at 96 GHz 
and the continuum flux density at 225 GHz are consistent 
with the 2 50 GHz flux density of S 250 GHz = 10.7 + 0.6 mJy 
JOmont et all2003l) . for a grey-body spectrum with a dust tem- 
perature of 45 K and a spectral index f3 = 1 .5. The correspond- 
ing far-infrared luminosity is Lfir ~ 4.3 x 10 13 L . The recent 
detection of J 1409+5628 at 350 yum is in agreement with these 
conclusions (Beelen et al. in preparation). 

3.2. Non-thermal radio emission 

Figures [5] and |4] show the VLA and VLBA radio images 
of J1409+5628 at 1.4 GHz. The VLA image at l'/6 reso- 
lution shows an unresolved source with a deconvolved size 
< 0'.'5, consistent with the upper limit derived from the CO 
measurements. The flux density at 1.4 GHz is Si.4ghz = 

0. 93 ± 0.022 mJy. The radio source position is 14 h 09 m 55 s .57 
+56°28'26'.'47 (J2000), and coincides within the astrometric 
errors (< 072) with the positions of the optical and the CO 
emission. 

The VLBA image at full resolution (12.5 mas x 6.7 mas, 
major axis position angle -9°; left panel of Fig.|4j has an rms 
noise of 46 fi]y beam -1 . There is no source brighter than Act in a 
field of 075 centered on the VLA position. We have also made 
tapered images at lower resolution and wider fields. At 34 x 
28 mas resolution the rms noise is 72 fjjy beam -1 . The tapered 
image is shown in the right panel of Fig. @] We still find no 
source brighter than 4<x within 073 of the field center. Near the 
center of the field there is a 'linear' feature extending about 07 1 
north-south with a peak surface brightness of 0.27 + 0.07 mJy 
beam -1 at 14 h 09 m 55 s :572, +56°28'26739 (J2000). It is possible 
that this linear feature corresponds to the J 1409+5628 radio 
source, but given the low surface brightness, we conclude that 
our VLBA observations can only set a Act upper limit to the 
surface brightness at 30 mas resolution of 0.29 mJy beam" 1 . 

4. Discussion 

4.1. The redshift of J1 409+5628 

The redshift derived from the CO measurements should 
correspond to the systemic redshift of the host galaxy of 
J 1409+5628, since it traces the extended molecular gas as- 
sociated with the quasar and not ionized gas undergoing en- 
ergetic processes linked with the AGN activity. The derived 
redshift for the molecular gas, zco = 2.5832 + 0.0001, is in 
the upper p art of the redshift range determined from the opti - 
cal spectra jKorista et alJll993UBarlow & Junkkarinenl[T9 94). 

1. e. 2.550 < z < 2.58. zco is close to the redshift of the 
stron g Ly-q absorption (z = 2.5758) seen in the Keck spec- 
trum JKorista etalJll993b . This absorption line is associated 
with narrow C iv absorption at 4352 A, indicating that the cor- 
responding gas is not part of the Lyman-Q' forest, but rather part 
of the host-galaxy. This is in agreement with the fact that this 
ioni zation state is not a ssociated with the central AGN (see, 
e.g.. lPetitiean et alJll994h and that its redshift is close to zco- 
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Fig. 3. VLA image of J 1409+5628 at 1.4 GHz at a resolution 
of 177 x 175 (FWHM), major axis position angle of 5°. The 
contours are a geometric progression in the square root of two, 
with the first contour level at 0.05 mJy/beam, and the highest 
one at 128 cr. Negative contours are dashed. 

The (systemic) velocity corresponding to zco is redshifted by 
1800 km s -1 from the z — 2.5624 high ionization ultraviolet 
emission lines such as C iv dKoristaetalJll993k a difference 
which has been found in many other quasars (Richards et alJ 
2002). 



lumi- 



4.2. Molecular gas emission 



The 



^CO(7- 



6) luminosity is about 1/3 of the L' CO(3 ^ 2) 
nosity (Table 0, a ratio which is comparable within a fac- 



>2 luminosity ratios measured in the 

IT 



tor of 2 to the 7-^6/3- 

high- z sources SMM J1401 1+0252 C0.2. iDownes & Solomon] 
2003), H1413+117, i.e. the C loverleaf (0.8 iBarvainis et alJ 
19971) . and J1148+5251 C0.6. iBertoldi et alJ l2003bl) . As in 
the other sources, the detection of CO(7— >6) emission in 
J 1409+5628 indicates the presence of warm and dense molec- 
ular gas. 

Assuming a constant brightness temperature from 
CO(3— »2) to CO(l— >0) and a conversion factor from L' CQ (1 ^ 0) 
to molecular mass of Xco ~ 0.8 M n (K km s" 1 pc 2 )' 1 
appro priate for ultraluminous galaxies ( Dow nes & Solo mon 
1998), the CO (3— »2) luminosity implies a mass of molecular 
gas of M(H 2 ) ~ 6.6 x 10 1() M Q . J1409+5628 appears therefore 
to be a system which is rich in molecular gas and comparable 
to the other high-z quasars detected in CO. 

4.3. Radio continuum emission 

Assuming a low frequency spectral index of -0.8, the 1.4 GHz 
continuum flux density implies a rest frame luminosity den- 
sity at 1.4 GHz of 4.0 x 10 25 WHz -1 . Using the far-infrared 
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Fig. 4. Left panel VLBA 1.4 GHz image at full resolution (12.5mas x 6.7 mas, major axis position angle -9°). The contours levels 
are: -0.3, -0.2, -0.1, 0.1, 0.2, 0.3, 0.4 mJy beam" 1 . Right panel The VLBA 1.4 GHz image tapered to 34 x 28 mas resolution. 
The contours are: -0.33, -0.22, -0.11,0.11, 0.22, 0.33, 0.44 mJy beam 1 . The central position corresponds to 14 h 09 m 55 s .5739, 
+56°28'26'/475 (J2000). 



luminosity definition of lCondonlll 19921), ex trapolated from the 
250 GHz flux density of Omont et al. (2003), a dust tempera- 
ture of 45 K and a spectral index ft — 1 . 5, the radio-to -FIR 
ratio is then q = 2.01, where q is defined bv lCondonl Jl992l) . A 
tight and linear correlation between radio and far-infrared lu- 
minosity has been found in star forming galaxies, w ith s mall 
scatte r over a few orders of magnitude in luminosity JCondonl 
Il992l) . The most recent consideration of this correlation found 
a value of q — 2.3 ± 0.3 for star fo rming galaxies from the 
IRAS 2 Jy sample JYun et alJboOll) . Hence, the value of q 
for J 1409+5628 falls well within the envelope defined by star 
forming galaxies. 



Cond on et al. I dl99lh derived an empirical upper limit to the 
brightness temperature for nuclear starbursts of order 10 5 K at 
8 GHz, while typical radio loud AGN have brightness temper- 
atures two or more orders of magnitude larger than this value. 
They also present a possible physical model for this limit in- 
volving a mixed non-thermal and thermal radio emitting (and 
absorbing) plasma, constrained by the radio-to-FIR correlation 
for star forming galaxies. For J 1409+5628, we use the observed 
surface brightness limit from the VLBA together with the mea- 
sured spectral index and redshift to derive a 4<x upper limit to 
the intrinsic brightness temperature of 2 x 10 5 K at 8 GHz typ- 
ical of star burst galaxies. This is in contrast with the results 
obtained by Momi ian et alJ J2004I) for a sample of three high- 
z quasars which were imaged with the VLBA. In these z > 4 
quasars, the radio-loud AGN dominates the radio emission on 
very compact size, i.e. a few milliarcsec, with intrinsic bright- 
ness temperatures in excess of 10 9 K. 



4.4. Starburst activity 

Overall, the properties of the CO emission and the radio emis- 
sion from J 1409+5628 are consistent with those expected for 
a star forming galaxy. The radio emission indicates a q value 
which is consistent with the radio-to-FIR correlation of star 
forming galaxies, and an intrinsic brightness temperature be- 
low 2 x 10 5 K at 8 GHz typical of nuclear starbursts. 

The ratio between the far-infrared and CO(l— >0) lumi- 
nosities of J 1409+5628 implied by the present observations, 
L FIR /Z^ O(1 ^ 0) * 500 L (K km s" 1 pc 2 )" 1 or L FIR /L C O(i^0) ~ 
1.1 x 10 7 is at the upper end of t he values derived for local ul- 
traluminous infrared galaxies bv lSolomon et alJ {l997), and in 
agreement with other high-z quasars where CO line emission is 
detected JCox et all2002j). 

Following Omont et the star formation rate 

(SFR) in J 1409+5628 can be derived from its far-infrared 
luminosity, assuming that star formation is the main contri- 
bution to the dust heating. Depending on the stellar initial 
mass function and the starburst age and duration, the ratio 
of SFR to infrared luminosity is in the range 0.8 to 2 x 
10 -10 M yr" '/Lq, yielding for J1409+5628 a SFR in between 
3 and 8 x 10 3 M yr _1 . The massive reservoir of molecular gas 
of 6 x 10 10 M can therefore fuel star formation at the rate im- 
plied by the far-infrared luminosity for about 10 million years, 
a typical duration for a starburst in local galaxies. The dense 
warm molecular gas is thus rapidly consumed unless its mass 
only represents a fraction of the total gas available in the host 
galax y of J1409+56 28 , perhaps in a halo, like in the case of 
M 82 ( W alter etalJ2002l) . Infalling gas could then further sus- 
tain the starburst activity in J 1409+5628 over longer periods 
of time. 
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Multiple structures such as double source, arc or ring, 
which are expected for gravitational lensing, are not detected 
in the maps of the non-thermal radio continuum on any scale 
from 10's of mas to a few arcsec. From the radio maps, there 
is therefore no evidence for lensing in J1409+5628. No correc- 
tion for magnification has therefore been applied in this paper. 
Deep optical and/or near-infrared observations would be useful 
to further search for evidence of lensing in J1409+5628. 

The VLBA observations, which resolve out structures 
larger than 0715, indicate that the 1.4 GHz radio continuum 
emission should be distributed on scales larger than 1 kpc. On 
the other hand, as shown by the CO observations, the molecu- 
lar gas is distributed in a region smaller than 5 kpc. If the radio 
continuum and CO emission are co-spatial, the present limits 
could be compatible with a scenario wherein the starburst ac- 
tivity in J1409+5628 is confined within a molecular torus or 
disk of size between 1 and 5 kpc. Similar structures with sizes 
of typically a few k pc are found in the ca ses of PSS 2322+1944 
JCarilli et alJl2003h and the Cloverleaf JVenturini & Solomonl 
120031) . 

Using the limits to the size derived from the CO and ra- 
dio observations and assuming that the molecular gas is dis- 
tributed in a disk, the dynamical mass is (3 < M^ ym < 
10) x \0 m sin~ 2 i M G , where i is the inclination to the line of 
sight. This range is comparable to the dynamical masses de- 
rived for other hig h-z quasars an d galaxies and to that of local 
m*-galaxies dGenzel et all2003l) . The ratio of molecular to dy- 
namical mass is therefore estimated to be (2.4 > MuJMd yn > 
0.6) sin 2 2, indicating that a large fraction of the total mass is in 
the form of molecular gas. 

5. Conclusions 

The millimetre and radio observations described in this pa- 
per provide strong evidence that the far-infrared luminosity of 
J 1409+5628 is related to massive star formation. First, the 
huge reservoir of dense and warm molecular gas (6 x 10 10 M ) 
can sustain the star formation rate of a few thousand solar 
masses per year implied by the far-infrared luminosity for 
about 10 million years. Second, the radio continuum emission 
is found to be consistent with the radio/far-IR correlation of 
star forming galaxies. The limits on the size of the CO and ra- 
dio continuum emitting region suggest that in J 1409+5628 star 
formation is taking place in an extended region, which could be 
a torus or disk-like structure not larger than 5 kpc. Observations 
of the CO emission at higher spatial resolution would be use- 
ful to further constrain the size of the star-forming region in 
J 1409+5628 and thereby its dynamical mass. 

The submillimeter and radio properties of the z = 2.58 
quasar J 1409+5628, which are similar to those seen in other 
high-z optically luminous quasars, are consistent with a coeval 
starburst and AGN. However, the relative importance of the 
warm dust heated by the star formation and the dust directly 
heated by the AGN is still poorly constrained. Observations at 
higher frequencies are needed to analyze the rest-frame near- 
to mid-infrared spectral energy distribution of distant quasars 
like J1409+5628. Such observations will allow us to derive the 
total infrared luminosity of high-z quasars, to estimate the con- 



tribution of the hot dust heated by the AGN, to better constrain 
the star formation history in their associated host galaxies and, 
thereby, to investigate further the coeval growth of stars and 
massive black holes and its evolution over cosmic time. 
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